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Abstract., - A new vehicle following controller . 
is proposed for autonomous intelligent vehicles. 

The proposed vehicle following controller not only 
provides smooth transient maneuver for unavoid- f 
able nonzero initial conditions but also guaran- . 
tees the asymptotic platoon stability without the 
- availability of feedforward information. Further- ^ 
more, the achieved asymptotic platoon stability 
is shown to be robust to sensor delays and an up- _ 
per bound for the allowable sensor delays is also : 
provided in this paper. ; 


I. Introduction 

Designing autonomous intelligent vehicles is important 
in the research of Advanced Vehicle Control Systems 
(AVCS) which is a major initiate in Intelligent Vehicle 
Highway Systems (IVHS). The main advantage of an au- 
tonomous intelligent vehicle is that it is considered as a 
“self-contained” system, i.e., it can operate together with 
other manually controlled vehicles without further techni- 
cal assistance from highway infrastructure. Since future 
Automated Highway Systems (AHS) is planned to evolve 
from today’s highway operation, the deployment of au- 
tonomous intelligent vehicles is of particular importance. 

An autonomous intelligent vehicle is assumed to be ca- 
pable of measuring (or estimating) necessary dynamical 
information from the immediate front vehicle by its on 
board sensors. The computer in the vehicle will then pro- 
cess these measured data and generate proper throttling 
and braking actions for controlling the vehicle’s move- 
ment. These longitudinal maneuvers must be performed 
as swiftly as possible within the rider’s comfort and safety 
constraints. 

Traditionally, vehicle following controllers are designed 
for single-mass (triple integration) models which do not 
account for any propulsion system dynamics, see, e.g., 
[1, 6]. In [8], Shladover included a simple first order en- 
gine model in the system dynamics and designed a linear 


vehicle following controller. It was shown that asymptotic 
platoon stability can be achieved by this linear controller 
when the drag forces (aerodynamic force and mechani- 
cal force) are neglected and the feedforward information 
is available. Based on the same vehicle model [8] with 
(nonlinear) drag forces taken into account, a nonlinear ve- 
hicle following controller was designed by Sheikholeslam 
and Desoer [7] using feedback linearization technique. In 
this case, asymptotic stability can also be achieved if the 
feedforward information is available. In [3], based on a 
more complicated vehicle engine model proposed in [5], 
Hedrick et al. proposed a sliding mode nonlinear con- 
troller to achieve vehicle following. The simulation results 
indicated that the controller has the potential of achiev- 
ing asymptotic platoon stability if the feedforward infor- 
mation is available. This observation was later verified 
with proof in [9]. In [4], Ioannou and Chien modified the 
nonlinear vehicle following controller proposed in [7] and 
showed that asymptotic platoon stability can be achieved 
by this modified controller without any feedforward infor- 
mation. This result enhances the feasibility of the future 
deployment of autonomous intelligent vehicles. 

In this paper, we propose a new vehicle following con- 
troller based on the nonlinear model proposed in [5] and 
[3] . The proposed vehicle following controller not only pro- 
vides smooth transient maneuver for unavoidable nonzero 
initial conditions but also guarantees the asymptotic pla- 
toon stability without the availability of feedforward infor- 
mation. Furthermore, we show that the achieved asymp- 
totic platoon stability is robust to sensor delay and an 
upper bound for the allowable sensor delays is provided. 

This paper is organized as follows. In Section 2 and 3 
, a vehicle longitudinal model and a safety distance pol- 
icy are briefly reviewed. In Section 4, we present control 
methodologies for two classes of nonlinear control systems 
based on the ideas developed in backstepping control tech- 
nique. Applying theses methodologies, we design vehicle 
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following throttle and brake controller in Section 5. The 
issues of designing asymptotic platoon stability and its 
robustness to sensor delays are discussed in Section 6. In 
Section 7, we use simulation results to demonstrate the 
effectiveness of our approach. At last, Section 8 gives a 
brief conclusion and possible future research directions. 

II. Vehicle dynamics model 

In this section, we introduce a longitudinal powertrain 
model for control system design. The derivation of the 
system dynamic equations is based on the following as- 
sumptions [9]: 

• Ideal gas law holds in the intake manifold. 

• Temperature of the intake manifold does not 
change. 

• There are no time delays in generating the power in 
the engine. 

• The drive axle is sufficiently rigid. 

• The torque converter is locked. 

• The brakes follow first order dynamics. 

The dynamics of the flow of air into and out of the intake 
manifold is described by 

7JT a — 771 a , — 771 a O 

where m 0 is the mass of air in the intake manifold and 
rn a i> ™<xo are the mass flow rates through the throttle valve 
and into the cylinders, respectively. 

Empirical equations developed for these flow rates are 

TTlai — rn a xPRi(m a )Tc{oc) 

771 a o — 771 ao ( > 771 a ) 

where m ax is a constant determined by the size of the in- 
take manifold; Tc() is the throttle characteristic, a non- 
linear function of the throttle angle a; Pri( ) is the pres- 
sure influence function describing the choked flow relation- 
ship. Notice that m a o is generally measured by steady- 
state engine tests and supplied in tabular form as a func- 
tion of the mass of air m a in the intake manifold and the 
engine speed w e . 

The engine’s rotational dynamics is given by 

Iw e = T n *t(u) c , m a ) - RT br - RrFtr ( 1 ) 

where I is the rotary inertia of the engine and the wheels 
referred to the engine side; R is the effective gear ratio 
from the wheel to the engine; Ti> r is the brake torque; 
T net is the net-engine torque which is also measured by 
steady-state engine tests and supplied in tabular form as 
a function of m a and the engine speed w e ; r is the effective 
tire radius; and F ir is the tractive force. 


The tractive force can be expressed as 
Ft r = K r sat ( i/J) 

where K r is the longitudinal tire stiffness; i is a constant 
determined by the road and tire condition (usually around 
0.15 [10]); sat(-) is the standard saturation function; and 
i is the slip between the wheels and ground given by 

V 

l = 1 - — — 

In addition, we adopt a linear brake actuator model 

Tbc ~ Pbr 
Tb 

where t\, is the actuator time constant, 7& r is the brake 
torque applied to the driven wheel and T& c is the com- 
manded brake torque. 

Finally, the longitudinal equation for the vehicle veloc- 
ity is given by 

Mv = Ftr — cv 2 — p>Mg (2) 

where cv 2 is the aerodynamic drag, fiMg is the rolling 
resistance, and M is the effective mass of the vehicle. 

Under the “no-slip” condition [9], i.e., 

v — Rrw e , 

equations (1) and (2) yield 

Jw c = Tnetfu/e, 77i a ) — cR 3 r 3 w 2 — RT br — 4>1 

where <pi = RrfiMg ; J = / - f Mr 2 is the effective inertia 
of the vehicle referred to the engine. 

With above discussions, the i tfl following vehicle has the 
following longitudinal dynamics, 


±i - 

v{ — Rrwe 

(3) 

W e — 

j[ Tnet(w c ,m a ) ~ cR 3 r 3 wl - RT br - fa) 

(4) 

TTlci — 

-77iao(w c ,m Q ) -1- m a xPRi(rn a )Tc{ct ) 

(5) 

II 

Tbc ” Tbr 

(6) 


Tf> 


where and Vi denote the position and velocity along 
the longitudinal direction. 

III. Safety distance policy 

For safe longitudinal operations, a following vehicle is re- 
quired to keep a safe distance from its preceding vehicle. 
From the traffic capacity point of view, the desired safe 
distance should be as small as possible. However, the ve- 
hicle’s performance capability, rider’s comfort constraint 
and other safety considerations impose minimum bound 
on this distance. In this paper, we will adopt a desired 
safety distance policy [4] for the ith following vehicle. 

Sd t = ^l( v ? — -}- A Vi + A3 (7) 

where Ai,A,A 3 are positive constants determined by the 
specified values of human reaction time, vehicle’s full ac- 
celeration and deceleration, and maximal allowable jerk 
during deceleration. 

While vehicle following is operating near a steady state, 
the velocity of the control vehicle is approximately equal 
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to the velocity of its preceding vehicle. Therefore, the 
safety distance policy can be well approximated by the 
constant time headway policy 

= Xvi + A 3 . (8) 

Let x» (x,_i resp.) and (i/,*_i resp.) be the position 
and velocity of the i ih ( i —I th resp.) vehicle. As shown in 
Figure (1), the spacing deviation for the ith vehicle from 
the desired safety distance is 

6{ := — 1 —xj — /j— Sd ^ 

= — t, — 1 1 — At;* — A 3 

where /,• is the length of controlled vehicle. 


Xi 

Vi 


Xi - 1 
Vi - 1 



<5* Sd, = Ai>j + A 3 


Figure 1: 

For a group of vehicles with each vehicle’s longitudinal 
dynamics described by (3) - (6), our control objective is to 
design a controller for each vehicle such that the following 
objectives are achieved: the spacing deviation Si can be 
regulated; the asymptotically platoon stability is achieved; 
and smooth transient response is guaranteed for non-zero 
initial spacing deviation and velocity deviation. 

To this end, it seems that input-output feedback lin- 
earization technique may provide a promising approach 
to deal with this nonlinear control problem based on 
the structure of the system. However, since the map- 
pings T ne< (-,-), and Pri( ) are supplied in tab- 

ular forms, their exact partial derivatives are not clearly 
identified. Consequently, feedback linearization method 
can not be applied directly. 

IV. Nonlinear Control Methodologies 

In this section, we will show how the basic ideas used 
in backstepping control design approach can be applied 
to controller design for two classes of nonlinear systems. 
The control methodologies developed will then be used to 
design vehicle following controllers in Section V.. 

Nonlinear control systems Class I 

Consider the following single-input single-output (SISO) 


nonlinear control system 


X = Jo(w) 

w = /l(u>t 2) 

Z = h(yu,z) + f 3^)91 (u) 

y = h(x,w y x m ,v m ) 


(10) 


where w,z £ R are state variables; y £ R is the output; 
/,' (i — 0,1, 2, 3) and g\ are smooth nonlinear functions; 
x m} v m £ R are bounded external signals; and u is the 
control input. 

The control objective is to design the input u so that 
the output y is regulated, i.e., lim^oo y(<) = 0, while 
the state variables z remain bounded. Our approach 
for finding an input u to achieved the control objective is 
based on the application of control Lyapunov function in 
the backstepping technique developed in [2]. 

The basic ideas of applying backstepping technique to 
the control design for system (10) are roughly summarized 
in the following. First, we neglect the dynamics of state z 
and treat z as the input, then find a control input z — Zd to 
achieve output regulation for the following reduced order 
system: 

x = f 0 (vu) 

w = ( 11 ) 

y = h(x, w, Vm) 

Second, construct a state feedback u from the computed 
Zd such that 

^lim (z(*) — Zd(t)) = 0 

Finally, we show the control objective is achieved for the 
closed loop system. 

For our approach, we make the following assumptions. 


Assumption 4.1 fs(z) is nonzero and ^f±(w,z) is 
bounded for all w,z £ R. 


Assumption 4.2 The system 

X = /i(r,u) 

y = x 

is bounded input bounded output (BIBO) stable 
where x, y, u E R. 


We now elaborate the control design procedure. Take 

V w (x,w,XmyVm) := h 2 (x, w t Xm , v m ), 7 i > 0 

as a Lyapunov function and evaluate the derivative of V w 
along the trajectory of (11). We get 


V w {x,w, X m ^ I'm) 


dh 

' ylh ~dt 


71 /i 

71 h 


dh . dh . 

W + o Xrn + O Vm 

OXm vVm 

dh , , . dh t . . dh . 

-r-fo{w) + — h (-w, z d ) + Xn 

OX OW OXrn 


\dh . dh . 
[d^ x+ d^ w 


dh . 

+ - — v m 
OVm 
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If, for k\ > 0, Zd is such that 

9 /i , dh . 

— ~Z x rn — — v m j 

OXm OVm 

then 

V w (x } W y Xm<,Vm) = “ ' 7i/ci/l 2 (x ) ti;, X m ,T; m ) 

For further developments, we will assume 

Assumption 4.3 There exists a z d satisfying (12) for 
all x, w , x m , v m in the domain of interest. 


Take 


Vu (x, Wy Zy Zd , Xm , %) 

:= V/, x m ,z;m) + “72 - z d) 2 > 72 > 0. 

as a Lyapunov function for (10). The derivative of V u 
along the trajectory of (10) is 


V u (x , W , Zy Zd y Xm , Vm ) 

dh dh . dh , 

+ “ Vm ] + 72 (^ z d){ z z d) 

OVm 

dh , , a/i , , . , dh . 

= ^ h lte Mw)+ dZ Mw ’* d)+ teZ Xm 

'm ] + -nh^-[fi(w,z) - fl( w ’ z d)\ 

dVm 

+72 (^ ” *d)[h( w i z ) + h{ z )9l( u ) *" *d] 
d h 

= -yikih 2 +'yih-—[h{w,z) - /i(^,2d)] 

3w 

+72(2 - Zd)[/2(w»«) + hi z )9l(u) - *d] 

If, for k 2 > 0, u is such that 

/3(*)si(v) = -k 2 (z - z d ) - h{v>yz) + z d 

- 7 7l^”^[/l(^.^) “ /l(^,2d)] 

72(2 - ^d) dw 


then 

Vu(x,u;, z, 2 <i,x m , v m ) — — — 72 ^ 2 ( z - z d Y 

Theorem 1 Consider the system (10) with the following 
proposed nonlinear state feedback controller 


u(Xy Wy Zy X m , V m ) 

= gf l (-j^j{—k2{ z — z d) — h(vjyz)-\-z d (12) 

where z d satisfies (12). Suppose that Assumptions 4.1 , 
4.2 and 4.3 are satisfied. Then for the closed loop system 
(10), (12), we have w remains bounded, y converges to 
zero and z converges to z d asymptotically. 


Proof: Let 


z z - z d 


Then the closed loop syatem (10), (12) yields a subsystem 

h — -fcl/l+ z d) “ Jl{™, z d)] ( 13 ) 

dw 

| — -k 2 z - — + Zd) - h{ w , z d)] ( 14 ) 

72 z dw 

By Assumption 4.1, we have 

Urn = *&(„, ,.) < 00 

Z-+Z d z — z d Oz 

This implies that ( h,z ) = (0,0) is an equilibrium of the 
system (13), (14). 

Take as a Lyapunov function for (13), (14). 

V(M) := h* + \ 72* 2 . 

which is a positive definite, descrescent, and radially un- 
bounded function. The derivative of V along the trajec- 
tory of (13), (14) is 

V — 7 i hh + 72 zz 

dh 

= + J\h — [fi(wyz) - fi(tUy z d ) ] - 72 * 2* 2 

dw 

-~nh^[fi(w,z) - h(vJyZ d )] 

dw 

— —7i fci h 2 — 72 k 2 z 2 

< 0 

Therefore, we see 

h, z £ L *2 n Loo * 

The boundedness of w can be established by the bounded- 
ness of z and Assumptions 4.2 and 4.3. Finally, from the 
well known lyapunov theorem, we conclude that h con- 
verges to zero and z converges to z d asymptotically. D 

Nonlinear control systems: Class II We now con- 
sider the nonlinear control system 

x = fo(w) 

w = fi(v)y z) + fi(n) 

z — f 2 (wyz) ( 15 ) 

h = f3(n) + 9i( u ) 

y — h(xyW,x m ,Vm) 

where x , w, zrj £ R are state variables; y £ R is the output; 
/, (i = 0, 1, 2, 3,4) and g\ are smooth nonlinear functions; 
x m) v m E R are bounded external signals; and u is the 
control input. 

The control objective is to design input u so that the 
output y is regulated while the state variables w, z, rj 
remain bounded. We assume 

Assumption 4.4 7^(77) is bounded. 

Assumption 4.5 The system 

X! = /i(xi,x 2 ) + fi(u) 
x 2 — f 2 (xi , X 2 ) 
y = [xi x 2 f 

is BIBO stable. 
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The control design for (15) is similar to the one for 
system (10). To start with, we neglect the dynamics of 
state 77 and treat 77 as the control input of the system (15), 
then try to find a control rjd to achieve control objective 
for the following reduced order system: 


X — 

foM 

(16) 

VJ = 

fl(W } Z) + A fad) 

(17) 

z = 

h(w t z) 

(18) 

y = 

h{x,Wy XmyVm) 

(19) 


Take 

V w (x,W f X m ,Vm) = ^pzh 2 (x,W y XmiVm), 73 > 0 

as a Lyapunov function and evaluate the derivative of V w 
along the trajectory of (19). We get 


V w (x, w, x m , v m ) 

, \dh . dh . dh , 

73/1 7-*+ + * x 

L ox ow oxm 


dh . 

t Vrn 

OV m 

= + + Mvd)] 

dh . dh . . 

+ a Xm + a Vm ■> 

OXjn OVm 


If, for > 0, Tjd is such that 


fi{Vd) = ~f\(vJ y z) + 
dh „ 




\ dxu J 
dh 


dx 


Io{vj) 


then 


■ o x ™ ~ t. — v m . 
OXm OVm 


V w (x, w, x m , v m ) = — rzhh 2 


( 20 ) 


Similarly, we assume 

Assumption 4.6 There exists an rjd satisfying (20) for 
all x , w, z , x m , i? m in the domain of interest. 

With Assumption 4.6, we take 

V u (x,W } r] i T} di Xm 1 Vm) = V w (x t W,Xm,Vm) + ^74 (*? ” 7?d) 2 

as a Lyapunov function and evaluate its derivative along 
the system (15). We have 

Vu(x, If/, 7J, TJd> Xm, Vm) 

= 73 M “/o(w) + |^[/l(v/,z) + A(7d) + Afr) - A fad)] 

dh , . 

+ - arm + Vm > + 74 (t? - T7 d )(7? - T} d ) 

OXm OVm 

dh 

= -73*3 h 2 + 73*— [A fa) - A fad)] 

ow 

+ 74(77 - 7?d)[/3(7?) + 91 fa) - T)d] 

If, for £4 > 0, u is such that 

5l fa) = ~*4 (77 - Vd) - A (T?) + Vd 

;4 rTS Aa(i) - /*(„„)] 

74(77 ~Vd) Ow 


then 


Vu(x, U/, 77, 77 d,X m , v m ) = ~ 73 * 3* 2 -74^4(77 - ? 7 d ) 2 


Theorem 2 Consider ike system (15) with the following 
proposed nonlinear state feedback controller 


u{x, W y Z,T7,X m , Vm) 

g( l ( -*4 fa - 7 u) - 73(77) + 77 j 

_J c73*f^[Afa) - / 4 (T7d)l ' 

74(77 - 7?d) OW 


( 21 ) 


where r]d satisfies (20). Suppose that Assumptions 4-4> 
4.5 and 4-0 are satisfied. Then for the closed loop system 
(15), (21), we have w,z remain bounded , y converges to 
zero and rj converges to 7}d asymptotically. 


Proof: The proof is similar to theorem 1. 


□ 


V. Vehicle Following Controller Design 

A vehicle following controller is required to maintain a 
desired spacing between vehicles and to guarantee asymp- 
totic platoon stability. The property that the spacing er- 
ror for a controlled vehicle can be regulated is referred to 
local stability. A platoon is asymptotically stable if there 
are no slinky-type effects [7] within a platoon. Researchers 
have found that local stability in vehicle following is not 
enough to guarantee asymptotic platoon stability. More- 
over, the unavoidable non-zero initial conditions occurring 
during various mode transitions, e.g., switching from man- 
ually control to automatic control, can generate transient 
torque large enough to degrade the driving quality. 

In this section, the control methodologies developed in 
Section IV.are applied to design a vehicle following con- 
troller with local stability and asymptotic platoon sta- 
bility. To deal with the undesirable transient response 
caused by non-zero initial conditions, we will filter the de- 
sired control effort by introducing an imaginary preceding 
vehicle in the controller design. Stability is guaranteed by 
the fact that the states of the imaginary preceding vehi- 
cle will converge to that of the true preceding vehicle ex- 
ponentially and the (imaginary) spacing deviation (from 
the desired spacing between the imaginary vehicle and the 
controlled vehicle) is regulated. With properly chosen de- 
sign parameters, the proposed controller achieves asymp- 
totic platoon stability which is robust to sensor delays. 

A. Controller Design 

The proposed controller is composed of a throttle con- 
troller, a brake controller, and a switching logic. The 
brake controller is to execute the decelerating operation. 
The throttle controller is to perform the accelerating and 
decelerating maneuvers while braking is not required for 
assistance. The switching logic is to properly activate and 
deactivate the throttle and brake controllers based on the 
needed control action at the current operating state. To 
be precise, the controller will continuously compute the re- 
quired throttle angle required by the control action. If the 
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calculated required throttle angle is greater than the min- 
imum throttle angle, say ao, the logic determines that the 
throttle controller alone is capable of handling the desired 
maneuver, and no brake torque is to be applied. If not, 
the logic will deactivate the throttle controller, i.e., keep 
the throttle angle at a 0 , and activate the brake controller 

to generate the proper brake torque. 

To smooth the transient response during vehicle ma- 
neuvering, we introduce for the i th (following) vehicle an 
imaginary preceding vehicle with dynamics characterized 
by the following equations 

Xi-l = Vi-l 

vi - 1 = -£ 2(^-1 - v»_i) - /h(x t _i - *»-i) / 22 x 

Xi— i (0) = (0) + li 4* Av, (0) -|- A3 

ui_l(0) = Vi(0) 

where i, t),_i can be viewed as the position and veloc- 
ity of the imaginary proceding vehicle for the i th vehicle; 

Pi = /?i(<5i(0), u*-i(0)— v»*(0)) and/? 2 = /? 2 (^(0), u t _i(0)- 
i>i(0)) are positive functions of <5,(0) and i(0) - v,-(0)) 
to be specified by designers. 

Remark 5.1 It is easily verified that if t>i_i = 0, i.e., 
the (true) preceding vehicle is traveling at constant 
velocity, it can be easily shown that (xj_i — Xj-iKO 
and converge to 0 exponentially. 

With suitably chosen parameters (3 1 and /? 2 , we can 
have proper convergence property of (x,_i — x,*_i)(f) 
and (£j_i - Vi-i)(t). 

Remark 5.2 Negative <5j(0) or v, — i (0) — t^i(0) may lead to 
the situation that the imaginary preceding vehicle 
is traveling ahead of the true preceding vehicle. For 
large negative value of (0) or Ui_i(0) — i>j(0), which 
is possibly an indication of impending collision, it is 
necessary to reflect this situation to the controller 
as soon as possible (which enables the controller of 
the controlled vehicle to be able to respond it prop- 
erly for avoiding collision). Therefore, the values of 
Pi and P 2 should be chosen in the sense that fast 
convergence rate is assured. 

Define 

Si := f,*_ 1 - Xi — Xvi - U - A3. (23) 

Compared (23) with (9), <5; can be regarded as the devia- 
tion of the desired spacing between the imaginary vehicle 
and the controlled vehicle. Furthermore, we see from (22) 

Si( 0) = 0. 

In order to shape the desired transient response, we adopt 
the idea of PID control and define a function to be regu- 
lated 

h := c p S x + Cj / Si d£ + (t>i_i - v<) (24) 

Jo 

where c p and cj are design parameters to be determined. 

The design of throttle and brake controllers are dis- 
cussed separately in the following. 


Vehicle following throttle controller 
Under the condition that the brake controller is deacti- 
vated, the vehicle longitudinal dynamic equations are re- 
duced to 


Xi = 

Vi = Rr We 

(25) 

W e = 

j[Tnct(we,Tn a ) - cRW e - <fc] 

(26) 

T 71 a — 

~rhao{w c , m a ) + m ax PRi(ma)Tc{c* ) 

, (27) 


We see that the system (25) - (27) with output function 
h given in (24) can be represented by equation (10) with 
the following variable and function substitutions 


(x,w, z } u) = (xij u»c, m a , O'), 
fo(w) = Rrw, 

fi(w,z) - ^j[Tnet(w,z) - cR 3 r 3 w 2 - 
h( w i z ) — -rh a o(w t z), 

h(z) = maxPRl(z), (x m ,v m ) = (x,_ 1 , 
t>*- 1), Si(u) = T c (u), 
h = (v m - Rrw) + c p (x m - x - \Rrw) 


+cj 


s: 


(x m - x - A Rrw)(£)d£ 


It is further verified that Assumptions 4.1, 4.2 are sat- 
isfied. Besides, the Assumptions 4.3 is also satisfied in 
the range of operation. By Theorem 1, we propose the 
following control law 


a 

marP^i (m,) {~^2( m a “ m ci,de$) + ?hao(^e i ™a) 
m-1/ . , *Yi (l + Ac*)Jlr N 

= T C ( +ma,d« + -f ) 

— [T ne t(w < , T7T a ) “ Tnet{vJe , ^T-a.dci]} 

m o 

where m a> des satisfies 

Tn aides ) 

= {A . - Vi) + ( c / + klCp)$i 

(1 + A Cp)Rr 

+k\cj I 6i d£ - P2 (vi— 1 - Vi) - Pi (r,_i - a:,) 

Jo 

+ki(v {- 1 - v^] + fa (28) 

From Theorem 1, it is clear to see 


Propositions Consider the system (25)-(27). The con- 
troller proposed in (28) - (28) will drive h to zero asymp- 
totically . 


While implementing the control law (28), m a>£ j ea is to be 
estimated by finite differencing sampling values of m a)des . 

We will delay the discussion of the convergence of 8i(t) 
until the brake controller is presented since in both control 
schemes we can show the same convergence property of 

m 


Vehicle following brake controller 
When the brake controller is activated, the throttle angle 
is kept at the minimum oro- In this case, the vehicle’s dy- 
namics is governed by equations (3) - (6) with a replaced 
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by the constant a 0 . Notice that the system (3) - (6) with 
output function (24) can be represented by (15) with the 
following variable and function substitutions 


(x,w,z y 77, u) = (xi t w Cl m a> T bry T bc ) y /o(^) = Rrw t 

1 R 

f\{w,z) - - [T ne t(w f z) - cR 3 r 3 w 2 - 4>i]> U(t}) = 

/2(u»,z) — — rh ao (w } z), (Xm , V m ) — (£«-l , V,_l) t 

fs{v) = - — »7i 31 M = “U, 

t 6 T b 

h = (v m — Hr-tt;) + Cp (r m — a: — Aftrty) 


+ c / 


/' 


(ar m — x — \Rrw)(£)d£ 


In addition, Assumptions 4.4, 4.5 are satisfied. And the 
Assumptions 4.6 is also satisfied in the operating range. 

To regulate the output function (24), we propose the 
following brake control law 


T bc 


T b[ ~Fi{T br ^br.des) 4" F br *f T br ^ e5 
Tb 

73 (1 + A c p )R 2 r 

h J 

74 J 


(29) 


where 


T^Pnetfwe, rn a ) - 4>l] — /*3 r ( 1 + x'e J> "'j [ c p(v»-l “ v ») 

+ (c 7 + kic p )$i + fcicj f* S t di - - Vi) 

-/? 2 (£i _ 1 - a;,) + - v;)] + 4>i 

By Theorem 3, we see 


From the definition of 6,- (9) and (30), we have 
(1 4- A Cp) 8 % 

= (1 + Ac p )(5;_i - vi - A Vi) 

~ (14- Ac p )v,_i - [-/? 2 (vi_i - vi- i) - Pi(vi-\ ~ ^i-l) 

+(cp 4- ki)(v t -i - v,) + (ki c p 4- cj)8i 4- kicj£] 
-A[-/? 2 (S,-1 “ V»-l) - - Vi-i) 

+(c p 4- fci)(v«_i — v») 4- (k\ c p 4- c/)6, 4- ^ic/5,] (31) 

Therefor, we have the following relationship: 

Mo 

[(i •+ Xc p 4 A0 a )» 3 + (ga 4 4 0iH t -i(O 

(1 + Ac,).* + (A C/ + Afcjcp + c p + fcx). 3 + (Afc lCj + c, + k lC p). + k lCj 
[A0 2 . 3 + (03 + A0i)* + 0l] w i-l(O 

^ (1 + Ac p )» 3 ■+ (Ac, + Xkicp + c p + ki)* 3 + (A*iC/ + c, + + kjc, 

( 32 ) 


Furthermore, from (22), we have stable transfer func- 
tion , v 

Vj-l(*) _ 023 + Pi / 33 ^ 

Vi — 1(«) S 2 +028 + 01 

From (32) and (33), we conclude that, by properly choos- 
ing design parameters c p , c/ and k 1 , we can make S{ con- 
verge to zero if v%-i is constant, (i.e., if the preceding 
vehicle is traveling at constant acceleration) and have sat- 
isfactory transient response of 

From the definitions of 8i (9) and 6* (23), we see 


Proposition 4 Consider the system (8)- (6) with a = c*o 
and output function ( 24 )• The controller proposed in (29)- 
(30) will drive h to zero asymptotically . 

Similarly, T& r ,des is to be computed numerically by finite 
difference sampling values of T brt des- 


Si = Si + Xi-x — Xi- 1 

As pointed out in Remark 4.1, £|_i — x»_i will converge 
to zero exponentially under the condition i*_i = 0. It 
follows that Si will converge to zero while the preceding 
vehicle is traveling at constant speed. 


Regulation of Si 

Recall that our goal is to regulate the spacing deviation 
Si in both throttle and brake control cases. This can be 
done by properly choosing control parameters c p) cj and 
k 1 as shown in the following. 

Let 

= *i. 

Since the engine/brake dynamics are much faster than 
the vehicle dynamics (which thus can be neglected in the 
stage of vehicle performance analysis), the vehicle dynam- 
ics of the closed loop system under either throttle control 
(28), (28) or brake control (29), (30) can be represented 
by 

vi — Rrw c 

_ _L — [-& 2 (vi-\ - vi- 1 ) - 0i (xi-i - x{- 1 ) 

1 4- Acp 

4-(cp 4- *i)(v»-i ~ v,) + (cj 4- k 1 c p )8 t 

+ k\cj f 6i dz] (30) 

JO 


VI. Asymptotic Platoon Stability 

In this section, we will show that by properly choosing de- 
sign parameters, the controller proposed in Section 5 can 
achieve asymptotic platoon stability when it is installed on 
each vehicle of a group of vehicles ( one following another) 
with safe distance rule (8). 

Asymptotic Platoon Stability 

Consider a group of vehicles all equipped with the pro- 
posed throttle controller (28) and brake controller (29). 
Since, at steady state of vehicle following, 

Vi - 1 = v,-_i and Si = <5,-, 
we see from (30) 
v* 

= + MK-i -vi)4-(c/4-*ic p )5i 

+* lC/ fo 6 '( z ) dz ] 

= c p 6 t 4- ( kicp + c/)5, 4- kicj J Q 6 t dt + fc 3 {vi-y - v,) 
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Differentiating equation (23) three times and substituting 
the derivative of by (34), we obtain 

"Si (0 

= Vi - 1 - Vi - X 

= c p 5i-i 4 (kjc p + c/)ii_i(t) 4 + ^i(v »_2 - vj-i) 

-[c p $, 4 (AlCp 4- cj)fc(t) 4 k 1 c I S % 4 fci(vi_i - v»)l 
— A[c p 6i +(k\Cp + cj)6 t + k\cjS t 4 *i(vi-i - #*)] 

From the above equation, we obtain the transfer function 
from 6i to <5j_ i 


(fcl 4 Cp)* 3 4 (fclgy 4 c/)» + fcjc; 

(1 + Ac p )i 3 +(*c/ + Xfcj^cp +c p + A!)* 3 + (XJc 1 ej -fcj + *1 <=*>)* +kic/ 

( 35 ) 

To avoid slinky-type effects, the disturbances caused by 
the lead vehicle in all frequencies should be attenuated 
along the following vehicles to insure that they do not be- 
come unreasonably large by the end. A sufficient condition 
for this to happen is for all i 


ci 



Figure 2: Parameter region for avoiding slinky effects 


| 5 'i , 0 | = |G 1 (ju;)|<l, for all tu > 0 (36) 

With G\(s) given in (35), the inequality in (36) yields 

[fcic; ~ (*1 + c p )u > 2 ] 2 + + cj ) 2 

I*l c / “ (* C J 4 + Cp + *i)ur 2 ] 2 + ui a [Afc 1 c/ + C/ - f kiCp) - (1 + Xc p )u< 3 ] 3 

< 1 for a.11 w > 0 

Simplifying the above inequality, we get 
(1 4 \cp) 2 w A 4 [A 2 c 2 4 A 2 * 2 c 2 
42 A/:i - 2(c/ 4 fcic p )]u/ 2 4 A 2 k^c 2 I > 0 

for all xv > 0 (37) 


Remark 6.1 When constant spacing safety policy (A = 
0) is adopted, inequality (37) for avoiding slinky- 
type effects reduces to 

v? — 2(cj 4 c p ) > 0 


Since c/ H- k\c p > 0 (to insure all the poles of 
Gi(s) are in the open left half complex plane), the 
above inequality can not be satisfied when w 2 < 
2 (c/ + fciCp). In other words, asymptotic stability 
can not be assured for low frequency disturbances 
under constant spacing safety distance policy. 


A sufficient condition such that (37) holds is 


A 2 c 2 4 A 2 kfc 2 p 4 2 A fciCp - 2 (c/ 4 *ic p ) > 0 
or equivalently 


where 


(<=/ - 4) 2 (Cp “ A a *?+«*l )2 _ 

A 2 + B 2 > 


2\ki 4 2 
A 4 (A*i 4 2) 
2Afci 4 2 
X 3 k\ (Xki 4 2) 


(38) 


Given k i > 0 and A > 0, the suitable values of parameters 
cj and c p satisfying inequality (38) reside outside shaded 
ellipse as shown in Figure 2. Consequently, if we choose 
ci,k\ y Cp outside the shaded ellipse as shown in Figure 2, 
asymptotic platoon stability can be assured. 


Asymptotic platoon stability under sensor delays 
In this subsection, the relationships between the sensor de- 
lays, the gains of the proposed controller, and the asymp- 
totic platoon stability will be investigated. The results 
obtained in this subsection can be used to quantify the 
performance requirements for the sensors for a specific de- 
signed controller. 

Let r be the time delay caused by the velocity sensor 
and the position sensor, such that the velocity and posi- 
tion terms in (28) and (30) are functions for t — r instead 
of t . Then the venicle dynamics of the closed loop system 
can be represented by 


v t (t) 


1 


1 4 c p x 
4*i cj 


[(c p 4 *i)(v;_i - v;)(* — r) 4 (ci 4 *i c p )S x (t - r) 
f 6i( (39) 

Jo 
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Differentiating both sides of(9) three times, we get 
‘Si (t) 

= \ Vi (t) 

= Cp 5*i - 1 (t - r) + (*ic p + c/)5*_i(< - r) 

-f*i cjdi — i (t — t) + *i — a — T ) “ vi—\ 0 — T )) 

— [c p 5j(i — r) -|- (kicp + c/)5i(t - r) 4- *ic/6,(f - r) 

+*i (v*_i (t — r) — t»i (t - t))] — A[c p 5 • (f “ T ) 

+ (*i Cp + cj)<5,(t — r) + *i c/5»(t — t) 

+*i(t5i-i(< - r) - Vj(t - r))] 

(40) 

Substituting (39) into (40) and taking Laplace trans- 
forms, we can derive the transfer function from 6j_i to 
Si 


Proof: Since the inequality in (44) can be rewritten as 

{[a 2 + 2a cos(iut) + l]w - 26 sin(u;r)}u; 3 
+{[6 2 - 2ac - 2c cos(u»r) - e 2 ]w + 2d s\n(wr)}w 
+(c 2 _ j2 _26d + 2de) > 0, 
asymptotic platoon stability is guaranteed if 

[a 2 + 2a cos(urr) + l]w - 26 sin(tuT) > 0, (46) 

[6 2 - 2ac - 2c cos(iar) - e 2 ]w + 2d sin(u/T) > 0, and (47) 
c 2 _ / 2 - 26d + 2 de > 0. ( 48 ) 

It is easily verified that 

(q-l) 2 _ (Ac p -i) 2 

26 2(A cj 4- A*j c p + c p + *i ) 

and 


(fcl + Cp)« 2 -KfclCy Cj)» + fcjCJ 

( e *r + \c p ).Z + (Xcj + Xfc lCp -I- Cp + fci)*2 + (Xfc icj + cj + kiCp). + fcic r 

(«) 

A sufficient condition for asymptotic stability is, for all i 

|G 2 (ju/)| 2 < 1, for allu/ > 0 

Substituting (41) into the above inequality, we obtain 

. (d-ew 2 ) 2 + f 2 w 2 | < 1? 

[d — bw 2 + (sinu/r)u ; 3 ] 2 + u > 2 [c — (a + cos(urr)u / 2 ] 2 
for allw > 0 ( 42 ) 

where 

a = Acp, 6 = Ac/ + A*i c p + c p + ki 

c = Afcic/ -f cj + *i c p , d = kicj ) (43) 

e = k\ 4" Cp, f — k\Cp cj 


fr 2 _ 2a c -2c- e 2 _ A 2 c 2 + A 2 fc?c 2 + 2Afc lC | - 2 (cj + fcic p ) 
2d _ 2*i cj 

such that condition (45) is equivalent to 


( ( a - l) 2 6 2 -2ac-2c- t 

< 2d 

Since a > 0 and 6 > 0, we see from (49) 

a 2 4- 2a cos(tur) + 1 > — 2 a + 1 > 26r > 26 

& 2 -- 2ac — 2c cos(tcr) — e 2 > 6 2 — 2ac — 2c — e 2 > 2dr 


(49) 


sin(rtt;) 

w 


> -2d 


sin(rti;) 


which guarantee the inequalities (46) and (47). Moreover, 


from (43), we see 

c 2 — f 2 — 26d + 2de = A 2 kfcj > 0 

which assures the inequality (48). 

VII. Simulation Results 


With equations in (43), condition for asymptotic platoon 
stability (42) is equivalent to 

[a 2 + 2a cos(tL>r) + l]u/ 4 - 26 sin(u/r)u/ 3 
+ [6 2 - 2ac - 2c cos(ujt) - e 7 ]w 2 + 2d sin(u/r)u/ 

+(c 2 _ f2 _ 2bd -f 2de) > 0, w > 0 (44) 

Proposition 5 Consider the vehicle longitudinal system 
(3) - (6) with control law (28), (29). The asymptotic pla- 
toon stability is guaranteed if 

. r (ACp ~ l ) 2 

r <C 2 (Ac/ + AA:ic p + c p + ki)' 

A 2 ^ + A 2 fe 2 c 2 + 2Afcic 2 - 2 (cj + fcicp) 

2kicj 

Remark 6.2 In Subsection 4.2, we have chosen \ 2 c] 4- 
\ 2 k\c 2 p + 2\h\c 2 - 2 (c/ + kic p ) to be positive to 
insure asymptotic platoon stability. Furthermore, 
Xcj 4 \k\c p + c p + ki and k x cj are also chosen to 
be positive to guarantee local stability (regulation 
of Si). Therefore, the right hand side of inequality 
(45) is positive. 


We consider vehicles following each other in a single lane 
with no passing. Each vehicle is assumed to be equipped 
with the proposed controller. The length of vehicles is 
assumed to be 4 meters. The following controller gains 
were selected for the simulations : 

Cp = 2, c/ = 0.5, c v = 2, 
k x = 5, fc 2 =40, *3 = 5, *4 = 1 » A = 1, A 3 = 2. 


Case Is Vehicle following with zero initial condi- 
tions: Six vehicles are assumed to follow each other and 

form a platoon in a single lane. The leading vehicle is as- 
sumed to accelerate from 9 m/sec to 15 m/sec, then to 21 
m/sec, and then to 27 m/sec. After achieves 27 m/sec, it 
then decelerates to 21.5 m/sec and then to 17 m/sec. Zero- 
initial conditions are assumed. The simulation results are 
shown in Figure (3). Good velocity tracking, small tran- 
sient spacing error and zero steady state spacing error are 
achieved for each vehicle. Moreover, no slinky- type ef- 
fects exist. In other words, asymptotic platoon stability 
is achieved. 
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Case 2 : Exit from the automatic lane: The follow- 

ing situation is considered : at time t = 0 sec, the leading 
vehicle changes lanes and the new vehicle target is 3.2 m/s 
faster and meters farther ahead than the previous one. In 
this situation, a suddenly change of the relative velocity 
and relative distance appears which is then confirmed by 
the on-board computer and the automatic control equip- 
ment is reset. Thus, non-zero initial conditions appear. 
The velocity, acceleration, and spacing deviation profiles 
shown in Figure (4) are quite smooth during the transient 
stage. 

VIII. Conclusion 

In this paper, we have studied the vehicle following con- 
trol problem for the autonomous intelligent vehicles under 
the constant time headway safety distance rule. Instead of 
using simplified linear vehicle following models frequently 
used in vehicle longitudinal control, we consider a nonlin- 
ear model that contains important attributes of engines 
dynamics. Using a newly developed nonlinear control 
technique, we are able to design throttle and brake con- 
trollers for the longitudinal control purpose with smooth 
maneuvers. One of features of this design is that the 
asymptotic platoon stability can be achieved with prop- 
erly chosen design parameters. We further show that this 
nice property is theoretically robust to a certain degree of 
sensor delays. The computer simulation results demon- 
strate the effectiveness of our control approach and en- 
hance the feasibility of practical AICC technology deploy- 
ment. 
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